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ABSTRACT 

In this work we show that Bardeen-Petterson accretion disks can exhibit unique, detectable features 
in relativistically broadened emission line profiles. Some of the unique characteristics include inverted 
line profiles with sharper red horns and softer blue horns and even profiles with more than 2 horns from 
a single rest-frame line. We demonstrate these points by constructing a series of synthetic line profiles 
using simple two-component disk models. We find that the resultant profiles are very sensitive to the two 
key parameters one would like to constrain, namely the Bardeen-Petterson transition radius rsp and 
the relative tilt j3 between the two disk components over a range of likely values [10 < rBp/(GM/c 2 ) < 
40; 15° < (3 < 45°]. We use our findings to show that some of the "extra" line features observed in 
the spectrum of the Seyfert-I galaxy MCG-6-30-15 may be attributable to a Bardeen-Petterson disk 
structure. Similarly, we apply our findings to two likely Bardeen-Petterson candidate Galactic black 
holes - GRO J1655-40 and XTE J1550-564. We provide synthetic line profiles of these systems using 
observationally constrained sets of parameters. Although we do not formally fit the data for any of these 
systems, we confirm that our synthetic spectra are consistent with current observations. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: 
Seyfert — line: profiles — X-rays: stars 



1. introduction 

Relativistically broadened emission lines have been ob- 
served in the X-ray spectra of over a dozen active galactic 
nuclei (AGN) (Mushotzky et al. 1995; Tanaka et al. 1995; 
Nandra et al. 1997). More recently similar features have 
also been discovered in Galactic black holes (GBHs) and 
black hole candidates (Miller et al. 2002a, b,c). These lines, 
from fluorescent Ka emission of iron, provide a unique 
diagnostic of the inner regions of accretion flows around 
black holes. Such fluorescent lines are produced when hard 
X-rays illuminate regions of optically thick disk material 
(George & Fabian 1991; Matt et al. 1991). Since the line 
energies are well known (6.4 keV for cold iron with ion- 
ization states below Fe XVII, 6.7 keV for Fe XXV, and 
6.9 keV for Fe XXVI) and the intrinsic widths are very 
small, the observed, broadened line profiles can provide 
direct information on the Doppler shifts and gravitational 
redshifts affecting the line-emitting material (Fabian et al. 
1989; Laor 1991). If one assumes the emitting gas to be in 
the form of a thin Kcplcrian disk, then these line profiles 
also provide information on the emissivity of the disk, the 
radial distribution of the emitting gas, and the observed 
inclination angle i of the disk (i — 0° is face on) (Laor 
1991; Pariev & Bromley 1998; Pariev et al. 2001). The 
profiles of these lines are also sensitive to features of more 
realistic disk models such as turbulent broadening and ra- 
dial inflow (Pariev & Bromley 1998; Armitage & Reynolds 
2003). 

The dependence of these line profiles on the radius of 
emission and inclination of the disk suggest that they may 
be useful as a diagnostic of warped or twisted accretion 
disks (Hartnoll & Blackman 2000). One such warping 



mechanism which is active at small radii around rapidly ro- 
tating black holes is caused by relativistic Lcnsc-Thirring 
precession. Here the frame dragging of the rotating black 
hole causes orbits that are inclined (or tilted) with respect 
to the spin axis of the black hole to precess. For test par- 
ticles, the precession frequency is SIlt — 2GMa/(c 2 r 3 ). 
For a tilted disk, the Lense-Thirring precession causes a 
radially dependent torque which, by itself, would cause 
the disk to twist up. However, in a disk this twisting 
must compete with whatever angular momentum trans- 
port mechanism (i.e. "viscosity") is driving the accretion. 
In such a case, the quasi-steady-state solution is character- 
ized by a two-component disk with the inner disk aligned 
with the symmetry plane of the black hole and the outer 
disk essentially retaining its original orientation (as set by 
the angular momentum of the gas reservoir). This is the 
so-called Bardeen-Petterson effect (Bardeen & Petterson 
1975). It is characterized in terms of a tilt angle (3 between 
the inner and outer disk components and a transition ra- 
dius tbp where the two components connect. 

If iron lines are emitted by the disk, then depending 
on where the lines are emitted relative to the Bardeen- 
Petterson transition radius, the resulting line profile may 
reveal unique detectable spectral features characterizing 
one or both disk components. First, if rpp lies within the 
line emitting region, part of the line flux will be emitted at 
one inclination relative to the observer while the remainder 
will be emitted at a different inclination. Since the blue 
wing (g max — Vmax/ve) of a broadened line profile is sensi- 
tive to the inclination of the emitting region, a tilted, two- 
component disk could produce multiple blue horns from 
a single (rest- frame) line. The red tail (g m i n = v m i n /v e ) 
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would also be sensitive to a Bardeen-Pctterson disk since 
it is responsive to the inner radius of the emitting region, 
which would obviously be different for an inner and outer 
disk component. The observer would see an integrated line 
profile that is in effect a superposition of the line profiles 
from the two disk components (minus any reduction in 
flux due to shadowing). Thus, as we will show, Bardeen- 
Pctterson disks can exhibit unique spectral features not 
seen in flat single-disk models, such as inverted line pro- 
files with sharper red horns and a softer blue horns or even 
profiles with more than 2 horns from a single (rest-frame) 
line. 

Perhaps the most important possible outcome of this 
research would be to directly constrain the scale of the 
Bardeen-Pctterson transition radius through observations. 
There is currently great uncertainty as to where the tran- 
sition radius is likely to lie, with estimates varying from 
rBP < 20r s (Nelson & Papaloizou 2000) to r B p > I00r g 
(Bardeen & Petterson 1975), where r g = GM/c 2 is the 
gravitational radius of the black hole. The expectation is 
that it will occur approximately where the rate of twist- 
ing by differential Lense-Thirring precession is balanced 
by the rate at which twists are diffused or propagated 
away by viscosity. The large uncertainty in r bp then is 
largely a product of our ignorance of the details of angu- 
lar momentum transport in accretion disks, particularly in 
tilted disks where the vertical shear component of viscos- 
ity may be crucial. Constraining tbp through observations 
would go a long way toward improving our understanding 
of black hole accretion disks. 

Even if rsp lies outside the line-emitting region, the 
presence of a Bardeen-Petterson disk could still be in- 
ferred in some cases. For instance, the tilted outer disk 
may shadow parts of the inner line-emitting region. The 
observable signatures of such shadowing was explored in 
Hartnoll & Blackman (2000), although for somewhat dif- 
ferent scales than are considered here. Furthermore, the 
Bardeen-Pctterson effect provides a natural explanation 
for why fits of some iron lines require inclinations that 
appear inconsistent with other observations of the same 
system. This is seen more often in GBHs where the incli- 
nation of the binary is often well constrained. For instance, 
the best-fit models for the iron line in V 464f Sgr give an 
inclination of i = 43° ± 15° (Miller et al. 2002b), whereas 
the optically determined inclination (of the far outer disk) 
is 60° < i < 71° (Orosz et al. 200f). Although these in- 
clination measurements are marginally consistent with one 
another, a true discrepancy of this type could easily be rec- 
onciled if one assumes that the iron line comes from the 
inner component of a two-component Bardeen-Petterson 
disk. 

We consider a range of values for the tilt (3 between 
the disks, with magnitudes in the range 0° < \(3\ < 45°. 
For GBHs, (3 < 30° is consistent with expectations based 
upon the formation avenues of compact binaries (cf. Frag- 
ile et al. 2001). ft is also consistent with the currently 
available observational constraints on such systems. For 
instance, as noted above, the best-fit models for the iron 
line in V 4641 Sgr give an inclination of i = 43° ± 15° 
(Miller et al. 2002b), whereas the optically determined in- 
clination is 60° < i < 71° (Orosz et al. 2001), yielding a 
possible tilt between the outer disk and the inner iron-line- 



emitting region of [3 > 2° (as we discuss below, the upper 
limit is uncertain because the orientation of the observer 
relative to the line-of- nodes is unknown). Interestingly, the 
radio jet observed from this source forms an angle with the 
line-of-sight of 6 < 6° (Orosz ct al. 2001). If the axis of 
this jet is coincident with the spin axis of the black hole 
or the angular momentum of the inner accretion flow, as 
is commonly assumed, then the tilt could be quite large, 
(3 > 54°. A similar, although less extreme, discrepancy 
is observed between the jet and disk inclinations in GRO 
J1655-40 (Fragile et al. 2001; Orosz et al. 2001), where 
the jet lies almost in the plane of the sky (9 = 85° ± 2°, 
Hjcllming & Rupcn 1995) and the disk has an inclination 
i '= 70.2° ± 1.9° (Greene et al. 2001), suggesting a pos- 
sible tilt of (3 > 15°. Another example is found in XTE 
J1550-564, for which the apparent supcrluminal motion of 
its radio jets (Hannikainen et al. 2001), constrains 6 to be 
< 53° while the observed inclination of the outer disk is 
i = 72.2°±5.2°, again suggesting a possible tilt of (3 > 14°. 

The case for Bardeen-Petterson disks around supermas- 
sive black holes is less clear. Nevertheless, many super- 
massive black holes are thought to be spinning rapidly 
(Elvis et al. 2002; Yu & Tremaine 2002) and two plausi- 
ble scenarios for supplying the necessary tilted accretion 
flow come to mind: 1) episodic accretion from tidally dis- 
rupted transiting objects, such as is thought to fuel the 
accretion of Sgr A*, the supermassive black hole at the 
center of the Milky Way; or 2) accretion following galac- 
tic merger events where the supermassive black hole has 
not had sufficient time to align with the accreting matter. 
Neither case would appear to have a preferred orientation 
and thus presents the possibility of a full range of tilts. 

The main goal of this paper is to highlight some of the 
unique features expected in relativistically broadened iron- 
lines emitted by Bardeen-Petterson disks. Although this 
is a highly idealized study, we show in §2 that the profiles 
are quite sensitive to the Bardeen-Petterson disk param- 
eters. Observations of such lines could therefore be used 
to set meaningful constraints on the physics of Bardeen- 
Pctterson disks. In §3, we proceed to discuss our results in 
terms of a few selected GBHs and AGN. We also illustrate 
that this technique for studying Bardeen-Petterson disks 
through relativistically broadened line profiles will be most 
successful when combined with other observations that can 
help constrain the system parameters. 

2. IMAGING AND LINE PROFILES OF 
BARDEEN-PETTERSON DISKS 

2.1. Disk Model 

We model the Bardeen-Petterson effect very simplisti- 
cally as two separate infinitesimally thin disks: an inner 
disk extending from r in to tbp and an outer disk ex- 
tending from tbp to r out . For most of our models we 
assume the black hole to be maximally rotating so that 
a/M — Jc/GM 2 = 1. In this case, the formal cut-off for 
the disk at the last stable orbit coincides with the hori- 
zon of the black hole, i.e. v ms 

= r B H = r„. However, we 
effectively truncate the disk at Ti n = l.05rsH by termi- 
nating photon trajectories that penetrate this radius. This 
is done for numerical convenience as the integration times 
of the photon trajectories become infinitely long as they 
approach the horizon. The choice of a/M = 1 was made 
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somewhat arbitrarily, and although this value cannot be 
reached by an astrophysical black hole, we demonstrate 
below that reducing a/M does not cause any more dra- 
matic change to our results than those attributed to other 
uncertainties in our medels. 

The outermost disk radius is set to r out = 100r g . The 
transition radius is varied over the range 10 < rsp/rg < 
40, consistent with the numerical results of Nelson & Pa- 
paloizou (2000). In such a simple two-component Bardeen- 
Pctterson disk model, we expect separate red and blue 
horns to be produced by each disk component. This is 
where the power of this technique lies. Since the blue horns 
will be sensitive to the relative tilt of the two disks and the 
red horns will be sensitive to the location of the transition 
radius, the full line profile gives information about the two 
parameters one would most like to constrain. 

All of the models in this work assume a twist-free 
disk (i.e. the linc-of-nodes of the two disks is a straight 
line). Based upon numerical simulations (Nelson & Pa- 
paloizou 2000) this seems consistent with expectations for 
a Bardeen-Petterson disk. We generalize our results by al- 
lowing the observer to be placed anywhere on a sphere of 
large radius (500r g ). The observer's location is described 
by two angles {9q, 4>q\, in the usual Boyer-Lindquist spher- 
ical coordinates centered on the black hole. Since the 
observer's location is often unconstrained, we consider a 
range of choices for # an d 4>o- However, as we show in 
§3.1, if the inclination of the inner or outer disk relative 
to the observer can be constrained, then the observer's 
location is no longer arbitrary. 

2.2. Ray Tracing 

To obtain an image of the accretion disk as seen by a 
distant observer, we integrate the geodesic equations of 
motion using a variant of the numerical ray tracing code 
described in Bromley et al. (1997). We follow the trajec- 
tories of photons from a pixelated grid of points at some 
large distance (r = 500r ff ) back to the accretion disk. The 
uniform grid may be associated with the field imaged by 
a CCD detector for instance. The observed frequency at 
each pixel is given by 



where subscripts o and e are observer and emitter, re- 
spectively, u is the 4-velocity of the emitter, and p is the 
emitted photon's 4-momentum. Generally, the emitter 4- 
velocity is specified by the disk model. In this work we 
simply use the Keplerian 4-velocity in the plane of the 
disk, ignoring effects such as radial transport and turbu- 
lent motion which have been explored in a general con- 
text elsewhere (Pariev & Bromley 1998). The photon 4- 
momentum is tracked while numerically tracing the pho- 
ton geodesic back in time from the observer's grid to the 
surface of the disk (inner or outer). We discard photons 
that arc intercepted by the black hole or otherwise miss 
the disk. Thus we properly accoint for the role of the 
black hole in absorbing photons emitted near the horizon. 
Here we assume that the transition region connecting the 
two disks is optically thin (i.e. photons passing between 
the two disk components near the transition radius don't 
"hit" anything). The validity of this assumption depends 
on the width of the transition zone and upon how much 



material on average is present within the transit region, 
neither of which is known at this time. 

In this work we also ignore the illumination of the disk 
by line photons emitted in other parts of the disk and the 
possible reflection of those photons to the observer. Be- 
cause of uncertainties in the emissivity function and the 
fact that the flux in the line is a small fraction of the X-ray 
continuum in this energy range, this is probably reason- 
able. Nevertheless, in future work we plan to explore the 
effect of line reprocessing in Bardeen-Petterson disks. 

2.3. Line Profiles 

From the pixelated disk image, we generate a synthetic 
line profile by making a histogram of the number of pix- 
els in different frequency bins. Each frequency bin is an 
accumulation of the contributions of individual pixels and 
thus represents an integral over the disk. The actual flux 
comes from weighting each bin by the emissivity e(r) and 
g 4 ; three factors of g coming from the relativistic invari- 
ant Ivjv^ ', where l v is the intensity, and the remaining 
factor arising because the line profile is an integrated flux. 
This weighting takes care of the effects of Dopplcr motions 
and gravitational redshifts. The emissivity e(r) is the local 
emissivity in the frame of the disk. Since some fraction 
of the photons emitted close to the black hole will be cap- 
tured by it, e(r) does not correspond directly to the energy 
received by the distant observer. 

Throughout this paper we cast our results in terms of 
a 6.4 keV rest-frame line energy. However, nothing in our 
methodology is tied to this value. In this sense our re- 
sults would be equally valid for any other spectral line 
that might be emitted from the region of the disk studied 
here. However, as will become clear, our results would be 
considerably confused if two or more lines were emitted 
with rest-frame energies separated by AE/E < 0.2. This 
could be particularly troublesome in accretion disks with 
a mixture of H-like, He-like, and cold iron. 

Throughout this work we assume an emissivity of the 
form e(r) oc r~ q and generally use g — 2.5. This is steeper 
than the r~ 2 expected for a centrally illuminating source, 
yet shallower than the r~ 3 expected from local energy dis- 
sipation within a disk. Obviously a shallower radial de- 
pendence would increase the flux contribution of the outer 
disk, whereas a steeper radial dependence would decrease 
it. Regardless, provided a detectable amount of flux is 
emitted by both disk components, our results will be quali- 
tatively insensitive to the choice of q since the locations (in 
frequency or energy) of the red and blue horns are insen- 
sitive to the specific flux distribution (Cunningham 1975). 
However, the strength of each feature, and hence its de- 
tectability, will depend upon the value of q. We explore 
this dependence further below. 

Other details of our results may also depend on the form 
of the emissivity function we have chosen. For instance, 
in the case of a central compact source (or in so-called 
"lamppost" models), the change in the angle of incidence 
of the illuminating radiation beyond tbp would cause a 
change in the slope of e(r). Accounting for such effects 
would require tracing the photon trajectories from their 
originating source, which is not done in the current work. 

Before considering specific black-hole systems, we first 
explore how our synthetic line profiles depend on the var- 
ious parameters of interest. We start by testing the de- 
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pendcncc of such profiles on the tilt angle (3 between the 
two disks. We consider four models with values of (3 in the 
range < [5 < 45°. For all of these models, the transi- 
tion radius is set to r BP = 15r g , the inclination between 
the observer and the outer disk is fixed at i ou t = 45°, and 
the observer is located at <fro = 90° (perpendicular to the 
line-of- nodes). Figure 1 shows both the image of the disk 
and the line profile for the (3 — 30° case. For illustration, 
we have plotted the line profiles of each disk component 
separately as well as the composit line profile. This gives 
a good feel for how the remainder of the line profiles are 
generated. Figure 2 shows the composit line profiles for 
the four models with varying (3. Clearly the /3 ^ cases 
are easily distinguishable. The most prominent signature 
of the Bardcen-Petterson profiles is the presence of an "ex- 
tra" horn, contributed by the inner disk component. Since 
the observer remains fixed relative to the outer disk in all 
these models, the blue edge g max , which is controlled here 
by the outer disk, remains fixed. However, the extra blue 
wing from the inner disk is sensitive to the tilt [3. In these 
models, a larger (3 translates into a smaller inclination for 
the inner disk which, in turn, results in the corresponding 
blue horn becoming redder. In fact, for the most tilted 
case considered, the blue horn of the inner disk (located 
at g « 0.84) is found at a lower energy than the red horn 
of the outer disk (located at g ~ 0.92). The key point, 
though, is that the separation of the various horns is a 
sensitive indicator of the tilt between the two disk compo- 
nents. In other words, measuring this separation provides 
a direct constraint on (3, although, as we show below, this 
conclusion does depend on the observer's location relative 
to the line-of-nodes. 

Next, we explore the dependence of a set of synthetic 
Bardeen-Petterson line profiles on the scale of the transi- 
tion radius rpp- For these models, (3 remains fixed at 30° 
and the observer is located at {0q = 15°, <f)o — 90°}, giv- 
ing inclinations of 15° and 45° between the observer and 
the inner and outer disk components, respectively. Fig- 
ure 3 shows the line profiles for five values of tbp in the 
range 10 < rsp/r g < 40. For these models, it is not the 
blue horn that shows the greatest variability, but rather 
the red horn. Thus, while the blue horn is sensitive to the 
tilt (3 between the disks, the red horn is sensitive to the 
transition radius rpp- This is because the red horn closely 
tracks the inner radius of a given disk component. Since 
increasing rpp is equivalent in our models to increasing 
the inner radius of the outer disk component, varying the 
transition radius causes a noticeable shift in the red horn 
of the outer disk from g w 0.76 for rpp = 10r g to g w 0.92 
for tbp = 40r 9 . Note that for the set of parameters con- 
sidered here, the inner disk does not exhibit a red horn; 
however, even if it did, it would not be very sensitive to a 
change in the transition radius since the minimum radius 
of the inner disk component would not change. Varying 
the transition radius also changes the amount of flux com- 
ing from the inner and outer disk components, so that the 
relative amplitudes of the various horns provides another 
diagnostic of Bardeen-Petterson disks. A final point to 
emphasize from Figure 3 is that the Bardeen-Petterson ef- 
fect causes detectable changes in the iron line profile even 
when rpp > 40r 9 . This is an important result, although 
as we show below, it is somewhat dependent on form of 



the emissivity function. 

Next we vary the black hole spin parameter a. Although 
the Bardeen-Petterson transition radius is expected to be 
a function of a, in this work we treat them as independent 
parameters so that tbp remains fixed at 15r g for the cur- 
rent set of models. The remaining parameters are also held 
at their canonical values ((3 = 30°, i ou t — 45°, </>o = 90°, 
and q = 2.5). Line profiles for four values of a/M between 
0.5 and 1 are shown in Figure 4. The biggest effect of 
changing a is that the inner radius of the disk ri n (assumed 
equal to r ms ) moves to larger radii for smaller a (r ms = 1, 
1.24, 2.32, and 4.23 for a/M = 1, 0.998, 0.9, and 0.5, re- 
spectively). Since photons are not emitted from as deep 
in the potential well for larger r^„ , the flux in the red tail 
is greatly diminished. Overall the flux contribution of the 
inner disk component is also reduced due to its decreased 
surface area. 

So far we have only considered observers with a line-of- 
sight perpendicular to the line-of-nodes (4>o = 90°). How- 
ever, the line profile for a given Bardeen-Petterson disk 
will depend on the specific location of the observer, so we 
now consider a set of line profiles made from the same 
Bardeen-Petterson disk configuration, but seen from dif- 
ferent viewing angles. To do this, we maintain a constant 
inclination between the observer and the outer disk while 
shifting the value of (j) . Effectively, the observer is cir- 
cling the angular momentum axis of the outer disk. For 
these models, the tilt, transition radius, and inclination 
of the outer disk relative to the observer remain fixed at 
(3 = 30°, tbp = 15r 9 , and i out = 45°. The observer is 
shifted through 180° from the positive t/-axis (0o = 90°) 
to the negative y-axis (<fio — —90°). Figure 5 shows the 
resulting line profiles. The most noticeable change is the 
dramatic shift in the blue horn attributed to the inner disk, 
which ranges from g w 0.81 for </> = 90° to g w 1.35 for 
0o = —90°. The contribution from the outer disk, on the 
other hand, remains roughly constant since the observer 
maintains a constant orientation with respect to this com- 
ponent and the transition radius remains fixed; therefore, 
the red and blue horns attributed to the outer disk are 
always located near g w 0.92 and g « 1.05, respectively. 
Nevertheless, there are still small changes in the outer disk 
contribution due mostly to shadowing by the inner disk. 
Thus far the results have shown that an iron line profile 
can be a sensitive probe for determining the relative tilt 
between two components of a Bardeen-Petterson disk and 
the radius at which the two components meet. However, 
these results have all assumed an emissivity of the form 
e(r) cx r~ q with q = 2.5. Although this is a reasonable 
choice based upon theoretical expectations and is consis- 
tent with some fitted line profiles, it is by no means the 
only plausible choice. Other values of q are certainly pos- 
sible, as are other functional forms for e (e.g. the "lamp- 
post" model is a popular alternative). We finish this sec- 
tion by considering a few other choices of q. Figure 6 
illustrates how the line profiles vary for choices of q from 
to 3. In all models, (3 — 30°, tbp = 15r 9 , i out — 45°, 
and 0o = 90°. For q < 2, the results are sensitive to our 
choice of r out , which remains fixed at 100r 9 in all cases. 
The profiles are normalized by the flux in the blue horn. 
It is clear from the figure that for q < 1 the large outer 
disk component swamps the emission from the inner disk. 
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Fig. 1. — Example disk image of a two-component Bardeen-Petterson disk with j3 = 30° and rgp = 15r 9 . The observer is located at 
: 90° (line-of sight perpendicular to the line-of-nodes) and is inclined 45° with respect to the outer disk and 15° with respect to the inner 



disk. 
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Fig. 2. — Synthetic line profiles for a sample of two-component Bardeen-Petterson disks with rgp = 15r 9 and various values for the tilt /3 
in the range < (3 < 45° . The observer is located at <j>0 = 90° and maintains a constant inclination of i ou t = 45° with respect to the outer 
disk. 
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Fig. 3. — Synthetic line profiles with the tilt fixed at (3 = 30° and the transition radius varied over the range 10 < rgp < 40. The observer 
is located at <j>o = 90° and is inclined 45° with respect to the outer disk and 15° with respect to the inner disk. 
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Fig. 4. — Synthetic line profiles with the tilt, transition radius, longitude of the observer, and inclination of the observer relative to the 
outer disk fixed at /3 = 30°, rgp = 15r 9 , <j>0 = 90° i ar >d iout = 45°, respectively. The emissivity assumes a form e(r) oc r~ q , with q = 2.5. 
The spin of the black hole is varied over the range 0.5 < a/M < 1. 
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Fig. 5. — Synthetic line profiles with the tilt, transition radius, and inclination of the observer relative to the outer disk fixed at /3 = 30° 



tbp = 15r s , and i ut = 45°, respectively. The longitude of the observer is varied over the range (a) < <f>Q < 90 and (6) —90 < 



< 0. 
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Fig. 6. — Synthetic line profiles with the tilt, transition radius, longitude of the observer, and inclination of the observer relative to the outer 



15r„ 



90°, and 



45° 



disk fixed at {3 = 30", Vbp 

value of q varied from to 3. All of these profiles are normalized using the peak flux in the blue horn at E 
are 2.73 X 10 4 , 50, 6.73, 1, and 0.157 for q = 0, 1.5, 2.0, 2.5, and 3.0, respectively. 
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Only for steeper power-laws does the inner disk contribu- 
tion become apparent. If we were to continue this exercise 
to consider even steeper power-laws, the inner disk compo- 
nent would eventually swamp the outer disk contribution. 
It is in the intermediate region, when both disk compo- 
nents are detectable, that this technique is most useful as 
a probe of Bardeen-Petterson disks. The exact range of q 
for which this will be true will depend on many parame- 
ters including r m i ni tbp, r maxi and the inclination of the 
observer relative to the two disk components. Neverthe- 
less, we have shown that this technique should work for a 
reasonable range of choices of all parameters. 

3. IMPLICATIONS FOR GALACTIC BLACK HOLES AND 
SEYFERT GALAXIES 

In the previous section, we showed that iron-line profiles 
could, in principle, be used as a diagnostic of Bardeen- 
Petterson disks. Useful science could be extracted from 
such profiles in the form of constraints on the relative tilt 
of Bardeen-Petterson disks and on the location of the tran- 
sition radius. However, the practical effect of including the 
Bardeen-Petterson effect in modeling of iron-line profiles is 
to add three new free parameters that must be fit: two new 
parameters come from the tilt f3 and transition radius tbp, 
themselves, and the third new parameter comes from re- 
placing a single inclination angle i with the two angles {9q, 
<f>o} necessary to describe the observer's location relative 
to the two disk components. Given the many uncertainties 
already inherent in fitting relativistically broadened lines, 
this is not an entirely attractive prospect. Nevertheless, 
we discuss some examples where there may be compelling 
reasons to make the effort. 

3.1. Galactic Black Holes 

Often the inclination i of an X-ray binary orbit can 
be determined through careful observation and fitting of 
the photometric light curve, the radial velocity curve, 
and other observables (e.g. rotational velocity of the sec- 
ondary), when the system is in X-ray quiescence. In cases 
where the primary black hole is a microquasar, the spin 
axis of the black hole can be inferred from the orientation 



of the relativistic jet, which makes an angle 9j et with re- 
spect to the line-of-sight of the observer. If both i and 
9 j e t are reasonably constrained, then the position of the 
observer (#0, </>o) m our models is confined to the family of 
points satisfying 

9q = 9j et 

cos i — cos P cos 9j et 



where we have utilized the law of cosines on a sphere. In 
this way, we can remove one of the extra degrees of free- 
dom required to fit a Bardeen-Petterson disk, leaving only 
the dependencies on (3 and tbp- At this time, this ap- 
proach probably gives the best hope for constraining the 
Bardeen-Petterson transition radius. (Spectral modeling 
of the continuum may also provide important constraints.) 
Following this procedure, we discuss two GBHs that are 
likely candidates for Bardeen-Petterson disks and provide 
sample synthetic line profiles of each. For now these syn- 
thetic line profiles are illustrative in nature since we do 
not fit actual observations. However, we do discuss the 
qualitative similarities between our synthetic lines and ac- 
tual line profiles in the literature (Balucinska-Church 2001; 
Miller et al. 2004). 

3.1.1. GRO J1655-40 

GRO J1655-40 is one of the best studied and most 
tightly constrained GBHs, with a binary period of 
2.62191 ± 0.00020 days, a primary (black-hole) mass of 
6.3 ± 0.5M Q , and a binary inclination of i = 70.2 ± 1.9° 
(Greene et al. 2001). It has also been previously identified 
as a likely candidate for a Bardeen-Petterson disk (Frag- 
ile et al. 2001). This characterization comes from noting 
the discrepancy between the binary inclination and the 
orientation of the highly collimated bipolar relativistic jet 
seen in this source during a 1994 outburst (9j et — 85 ± 2°; 
Hjcllming & Rupen 1995). If we assume that the orienta- 
tion of the jet is fixed by the angular momentum axis of the 
black hole, then the tilt between the binary orbital plane 
and the black hole symmetry plane is (3 > 10.9°. Due to 
the remaining degree of freedom to select the precise loca- 
tion of the observer (i.e. choosing <fio in our models), we 



Fig. 7. — Disk image of GRO J1655-40 for an observer at {9q, </>o}={85°, 32°} (/3 = 30°). Using observational constraints we fix the 
inclinations of the inner and outer disk components at jj n = 85° and i ut = 70.2°, respectively. The transition radius is set at rgp = 15r g . 



can not place a tight upper limit on the tilt. The formal 
upper limit is (3 < (180 - 10.9)° = 169.1°. Nevertheless, 
we argue on physical grounds that (3 < 30° is probably a 
reasonable upper limit. For instance, values of (3 > 90° 
would imply that the inner disk is orbiting in the opposite 
sense of the outer disk. Although such large tilts may be 
possible in some systems (e.g. systems that have under- 
gone binary capture or binary replacement), such counter- 
rotating disks seem a highly unlikely scenario. In fact, it 
is hard to imagine that a Bardeen-Petterson disk struc- 
ture could be supported if (3 were much larger than about 
45°, as this would require a dramatic reorientation of the 
angular momentum of the gas in the transition region. 
Adopting a limit of f3 < 30° and using equation 2 to fix 
the observer's location {8o, </>o}, we consider two possible 
tilt angles for GRO J1655-40: j3 = 15° and 30°. The corre- 
sponding observer locations are {#0, 4>o}—{85°, 81°} and 
{85°, 32°}. In both cases, the inclinations of the inner 
and outer disks are ii n — 85° and i out — 70.2°, respec- 
tively, consistent with the observations. Thus the inner 
disk is viewed more edge on than the outer disk, although 
both are viewed from fairly low latitudes, as shown in the 
disk image in Figure 7. Notice that the inner disk is viewed 
from such a low latitude that photons from the bottom of 
the disk are bent enough to propagate through the gap 
at the transition radius and make it to the observer. Al- 
though this is an interesting illustration of light-bending 
effects, it also highlights the limitations of our assumptions 
of an infinitesimally thin disk and an optically thin tran- 
sition radius. In a real system, these photons would likely 
be blocked from view by the flaring of the disk or optically 
thick gas in the transition region. Nevertheless, we include 



these photons in the synthetic line profiles shown in Fig- 
ure 8. The four lines in the figure correspond to the two 
observer locations and two values for the transition radius: 
tbp = {15, 30}r s . All of the models assume a maximally 
rotating black hole, an outer radius of r max — 100M, an 
observer at r = 500 M, and a disk emissivity of the form 
e oc r -2 - 5 . The noise in the profiles is a result of the poor 
resolution of the inner disk, which again is seen nearly 
edge-on. 
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The peak in the synthetic line profile at E w 7 keV is 
attributable to the blue horn of the outer accretion disk 
component overlaid on the much broader line from the in- 
ner disk. The energy of this peak is very close to the peak 
line energy seen in the ASCA observations of this source 
(Miller et al. 2004). Assuming a rest-frame line energy 
of 6.4 keV, this observed line energy requires a disk in- 
clination of i > 45°, consistent with all other indications 
that this system is seen at high inclination. Rest-frame 
line energies of 6.7 or 6.9 keV (from highly ionized iron) 
could be fit with smaller inclinations, but would require 
the emitting region to be far from the hole. 

The extremely high inclination of our model inner disk 
leads to a very large value for the blue maximum {g max = 
1.5 or E max = 9.4 keV for a 6.4 keV line), below which 
there is a shallow blue plateau. Although such high-energy 
line emission is not apparent in the observed spectrum of 
GRO J1655-40 (Miller et al. 2004), this does not necessar- 
ily rule out a highly inclined inner disk component. Such 
a component could be present, yet have its emission ob- 
scured through self-absorption by the finite-thickness of 
the disk itself (unlike the infinitcsimally thin disks con- 
sidered here). Or it could be that the high-energy line 
emission is hidden by absorption features. This is a po- 
tential problem for all relativistically broadened lines, not 
just those from Bardeen-Petterson disks. Absorption is 
particularly problematic at energies above about 7 keV. 

3.1.2. XTE J 1550-564 

XTE J1550-564 has many similarities to GRO J1655- 
40, although its parameters are generally not as well con- 
strained. Despite the larger uncertainties, it is still clear 
that the system hosts a black-hole primary (8.36M Q < 
Mi < 1O.76M ). The binary period (1.552 ± 0.010 days) 
and inclination (67.0° < i < 77.4°) are also reasonably 
well known. Similar to GRO J1655-40, XTE J1550-564 is 
a microquasar. During a 1998 outburst, large scale rela- 
tivistic jets were observed in both X-ray and radio wave- 
lengths (Hannikainen et al. 2001; Corbel et al. 2002). The 
radio components displayed apparent separation velocities 
of > 2c (Hannikainen et al. 2001). Although the exact 
orientation of the jets was not confirmed, the apparent 
supcrluminal motion itself constrains the angle of the jet 
relative to the linc-of-sight to be 6j et < 53°. Similar to 
GRO J1655-40 and V 4641 Sgr, this jet orientation is in 
disagreement with the inclination of the binary orbit; in 
this case, the implied tilt is (3 > 14°. Although the limit 
on this tilt is very similar to the limit for GRO J 1655-40 
and the binary inclinations are quite similar, there is a 
key difference between the two systems: for GRO J1655- 
40, the implied inner disk is viewed more edge-on than the 
outer disk, whereas for XTE J1550-564, the orientation is 
reversed and the outer disk is viewed more edge-on. Since 
the superluminal motion of the jet only gives us an up- 
per limit on 0, we are unable to use equation 2 to further 
constrain the observer location for XTE J1550-564. Nev- 
ertheless, for illustration, in Figure 9 we assume a jet angle 
9 j e t = 45° and a tilt (3 = 30°. The corresponding observer 
locations is {9o, </>o}={45°, —60°}. The inclinations of the 
inner and outer disks are ii n = 45° and i ou t = 72.2°, re- 
spectively. Again we assume a maximally rotating black 
hole, an outer radius of r max — 100M, an observer at 
r = 500M, and a disk emissivity of the form e cx r -2 ' 5 . 



We consider two values of tbp = {15,30}r g . Since the 
outer disk component for XTE J1550-564 is nearly identi- 
cal to what we used for GRO J1655-40, it is not surprising 
that its contribution to the line profiles is quite similar 
(compare Figures 10 and 8). The differences between the 
two figures are mainly attributable to the inner disk com- 
ponent and to differences in the shadowing of the various 
models. The synthetic spectrum of XTE J1550-564 in Fig- 
ure 10 has a peak at E « 7 keV. Some of the photons in 
this peak come from the outer edge of the inner disk com- 
ponent (ii n = 45°), while the rest come from the inner edge 
of the outer disk (i out — 72.2°). This value is roughly con- 
sistent with the peak of the observed line in XTE J 1550- 
564 at E w 6.8 keV (Miller et al. 2004). For a rest frame 
energy of 6.4 keV, this requires a disk inclination i > 45°. 
For tbp = 15r g , the synthetic spectrum also shows some 
excess emission in the range E = 7 — 8 keV. This comes 
from the inner edge of the outer disk component and is 
not seen for tbp > 30r g since the Keplerian velocity of 
the orbiting gas is not high enough in such a case. In- 
terestingly, excess emission above the continuum fit also 
appears over the same energy range in the observations 
of XTE J1550-564 (Miller et al. 2004). Line flux in this 
energy range requires a disk inclination of i > 60° , consis- 
tent with the binary inclination and the inclination of the 
outer disk component in our models. Although we have 
not formally fit the data for XTE J1550-564, it is clear that 
the observed profile is consistent with a Bardeen-Petterson 
disk configuration. 

3.2. Seyfert Galaxies 

Although GBHs may provide the best constrained sys- 
tems for studying iron-line profiles, at this time their data 
do not have as high of a signal-to-noise value as the obser- 
vations of Seyfert 1 galaxies. The high quality of Seyfert 
1 data makes it appropriate to at least consider this class 
of object as part of this study. Furthermore, there is a 
potential puzzle that a Bardeen-Petterson disk might help 
rectify: According to the unified model of active galactic 
nuclei, type 1 Seyfert galaxies are believed to be viewed 
at low inclination angles (i < 30°) in order that they not 
be obscured by the dust torus presumed to surround the 
accretion disks at large radii. However, such low inclina- 
tions would prohibit the blue wing of an iron line from 
extending much above its rest-frame energy (g max < 1). 
Therefore, the presence of detectable emission above 6.4 
keV is a potential problem. It has usually been interpreted 
as evidence that the emission is coming from highly ion- 
ized states of iron (Fe XXIV- XXVI). However, we consider 
a Bardeen-Petterson profile such as the ones shown in this 
work to be a potentially more conservative, alternative ex- 
planation. In the context of the unified model, the outer 
disk of a Seyfert 1 galaxy would be required to have a low 
inclination, but the Bardeen-Petterson model would allow 
the inner disk to take on a broader range of inclinations. 
In particular, a higher inclination for the inner disk would 
allow the blue wing to extend above 6.4 keV without re- 
sorting to higher ionization states. 

3.2.1. MCG-6-30-15 

One of the best studied iron-line sources is the nearby 
type 1 Seyfert galaxy MCG-6-30-15. It was the first 
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Fig. 8.— (a) Synthetic line profiles for GRO J1655-40 for an observer at {6> , </>o}={85°, 81°} (B = 15°) and (b) {6 , O }={85°, 32°} 
[P = 30°). In all cases, we assume the black hole is maximally rotating, r m ax = lOOAf, and e oc r~ 5 . In each case, we present profiles for 
two values of the transition radius (rsp = {15, 30}r 9 ). 
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Fig. 9. — Disk image of XTE J1550-564. We use observational constraints to fix the inclinations of the outer disk component at i ou t = 72.2°. 
We assume the inner disk is inclined ij n = 45°, consistent with the observational constraint i< n < 53°. The transition radius is fixed at 
tbp = 15r 9 and the observer is located at {9o, </>o}={45°, —60°}. 




Fig. 10. — Synthetic line profiles for XTE J1550-564. For these, we assume the black hole is maximally rotating, r max = 100M, ear 2o , 
and /3 = 30°. We present profiles for two values of the transition radius (tbp = {15,30}rg). 



12 



source to show a relativistically broadened iron-line fea- 
ture (Tanaka et al. 1995) and remains one of the most 
cited examples. In fact, the data for this source are of 
such high quality, they actually reveal a complex set of 
features that can not be adequately fit as a single line 
(Fabian et al. 2002; Ballantyne et al. 2003). The key fea- 
tures include a strong, broad peak near 6.4 keV, a very 
extended red tail down to E < 3 keV, and an "extra" 
horn or peak near 6.9 keV. As mentioned above, a 6.9 keV 
emission feature is a potential problem for Seyfert galax- 
ies in the context of the unified model since disks seen at 
low inclinations (like Seyfert- Is) can not produce strongly 
blue-shifted lines. Thus this feature is often attributed to 
highly ionized states of iron (Fe XXIV-XXVI). Alterna- 
tively we show that these complex features can be qualita- 
tively explained using a "single" 6.4 keV iron line from a 
Bardeen-Petterson disk. The key principle is that we can 
satisfy the unified model of AGN provided we model the 
outer disk component with a small inclination while we are 
free to choose a broader range of inclinations for the inner 
disk. We suggest therefore that, in the case of MCG-6- 
30-15, it is the outer component that produces the core 
line near 6.4 keV while an inner disk component produces 
both the extended red wing and the 6.9 keV feature. In or- 
der to do so from a 6.4 keV rest-frame line, the inner disk 
must be inclined at least 45° with respect to the observer 
(Pariev et al. 2001). If we assume that the outer disk is 
inclined 30°, this gives a tilt (3 > 15° depending on the 
angle between the line-of-sight and line-of-nodes ([3 = 15° 
if 0o = 90°; (3 > 15° if <t> ^ 90°). 

In Figure 1 1 we show a synthetic line profile that quali- 
tatively matches the data for MCG-6-30-15 (Fabian et al. 
2002). This profile assumes the observer's line-of-sight 
is perpendicular to the line-of-nodes (<j>o = 90°) so that 
the tilt is (3 — 15°. Unlike most of our previous mod- 
els, in which the black hole was maximally rotating, here 
we assume a spin a/M = 0.5, giving an inner disk ra- 
dius of r m i n = A.2M. The transition radius is taken to 
be tbp = 20M. As in previous models, we assume an 
emissivity of the form e oc r" 2 - 5 . We note, however, that 
this emissivity may not be steep enough to reproduce the 
observed red wing flux. Nevertheless, it is worth point out 
that our relatively simple physical model is successful in 
matching many of the unusual features of the MCG-6-30- 
15 line profile. As described earlier, a nice property of our 
model is that it provides a natural explanation for why the 
flux in the line core of MCG-6-30-15 does not match the 
flux in the red wing: these two flux contributions come 
from different disk components with different inclinations 
relative to the observer (in this case, the outer disk is seen 
more face on). 

4. CONCLUSIONS 

We set out in this work to demonstrate that relativis- 
tically broadened iron lines can reveal important details 
about Bardeen-Petterson disks in GBHs and AGN. In 
support of this goal, we showed that such spectral lines 
are sensitive to the two parameters of greatest interest 
in Bardeen-Petterson disks: the tilt between the two disk 
components and the transition radius where they meet. At 
this time, this approach probably gives the best hope for 
determining these parameters, although spectral modeling 



of the continuum can also provide some constraints. 

We applied our results to a small set of potentially inter- 
esting systems. One class of such systems is microquasars. 
These objects are often better constrained in terms of the 
relevant parameters than other potential targets such as 
AGN. Any fits to iron-line profiles in microquasars must be 
consistent with observations of the binary orbital parame- 
ters, the orientation of the relativistic jet, and spectral fit- 
ting of the thermal (accretion disk) component. This kind 
of consistency check will be particularly helpful in ruling 
out competing effects. Unfortunately, current observations 
of iron lines from microquasars do not have sufficiently 
high signal-to-noise to resolve the extra features that are 
characteristic of Bardeen-Petterson disks. Nevertheless, 
current or future missions should be able to achieve the 
required sensitivity. In anticipation of this, we have pro- 
vided predicted spectra of two likely Bardeen-Petterson 
microquasar candidates: GRO J1655-40 and XTE J1550- 
564. 

Although AGN are generally not as well constrained as 
microquasars, they nevertheless represent another impor- 
tant class of objects in the context of this work. Further- 
more, they have historically yielded higher quality line pro- 
files. In this work we focused on one example: MCG-6-30- 
15. Here we showed that a Bardeen-Petterson disk profile 
could qualitatively explain the complicated line profile ob- 
served in that system, without resorting to extra lines from 
other ionization states of iron. Although other models 
have been proposed that produce similar features, we feel 
the Bardeen-Petterson explanation is an intriguing possi- 
bility. The discovery of a Bardeen-Petterson disk in an 
AGN would reveal a wealth of information about the en- 
vironment of the accretion flow and the history of the sys- 
tem. 

The most restrictive limitation of this technique is that, 
to be fully exploited, both components of the Bardeen- 
Petterson disk must produce detectable levels of line flux. 
Provided this condition is met, we have shown that the line 
profiles exhibit unique signatures, such as inverted (red- 
horn-dominant) lines or lines with more than two horns, 
that are not easily mimicked by other effects. Nevertheless, 
there are competing effects one must be aware of. Prob- 
ably the most important are line blending, for instance 
from low and high ionization states of iron, and disk warp- 
ing due to mechanisms other than Lense-Thirring preces- 
sion. Line blending will be particularly confusing if two or 
more lines are emitted with rest-frame energies separated 
by AE/E < 0.2. In such cases, one will have to rely on 
other observational and physical constraints of the system 
to disentangle the possibilities. 

In closing, we would like to suggest that it might be 
fruitful in the future to conduct a systematic search for 
Bardeen-Petterson features among all observations of rela- 
tivistically broadened iron lines. Such a search could yield 
important statistical bounds on questions such as: How 
many systems have tilted accretion disks? What is the 
average tilt in such systems? What is the average transi- 
tion radius? Answering these questions will give important 
clues toward understanding accretion-disk physics and the 
evolutionary history of accreting black holes. 

P.C.F. gratefully acknowledges useful discussion and 
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Fig. 11. — Synthetic line profile providing a qualitative fit to MCG-6-30-15. For this model, the disk parameters are r m i n = 4.2M, 
rgp = 20M, fi = 15°, r max = 100M, and e tx r~ 2 ' 5 . The observer is located at ipo = 90° and has inclinations of i ou t = 30° and ii n = 45° 
with respect to the outer and inner disk components, respectively. 
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